Two-mode single-atom laser as a source of entangled light 



M. Kiffner/'B M. S. Zubairy, 1 ' 2 >Q J. Evers, 1 '! and C. H. Keitel 1 '! 

1 Max-Planck-Institut fur Kernphysik, Saupfercheckweg 1, 69117 Heidelberg, Germany 
2 Institute for Quantum Studies and Department of Physics, Texas A&M University, College Station, Texas 77843, USA 

A two-mode single-atom laser is considered, with the aim of generating entanglement in macro- 
scopic light. Two transitions in the four-level gain medium atom independently interact with the 
two cavity modes, while two other transitions are driven by control laser fields. Atomic relaxation 
as well as cavity losses are taken into account. We show that this system is a source of macroscopic 
entangled light over a wide range of control parameters and initial states of the cavity field. 
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PACS numbers: 03.67.Mn, 42.50.Dv, 42.50.Pq 

I. INTRODUCTION 

Quantum entanglement is known to be the key re- 
source in many applications of quantum information and 
quantum computing These phenomena range from 
quantum teleportation [2, S| and quantum cryptography 
4] to quantum implementation of Shor's algorithm [5] 
and quantum search [6(. It is therefore not surprising 
that there has been a great deal of interest in the gener- 
ation and measurement of entanglement in recent years. 

Entangled states have been considered traditionally 
between individual qubits. However, it has been shown 
that continuous variable entanglement can offer an ad- 
vantage in some situations in quantum information sci- 
ence [jj. One reason for this is that continuous vari- 
able entanglement often can be prepared uncondition- 
ally, whereas the preparation of discrete entanglement 
usually relies on an event selection via coincidence mea- 
surements. The classic scheme for the generation of con- 
tinuous variable entanglement is the parametric down- 
conversion. Starting with the first demonstration by Ou 
et al. {&], the generation of entanglement in such sys- 
tems has been achieved in several experiments Q- It 
still remains, however, a challenge to generate entan- 
glement in macroscopic light rather than on the few 
photon level. Promising candidates for the generation 
of macroscopic lig ht entanglement are optical ampli- 
fiers [3, E3, EH, E2 [HI]. For example, it was shown re- 
cently that a two-mode correlated spontaneous emission 
laser (CEL) [II El 

can lead to two-mode entanglement 
even when the average photon number in the field modes 
are very large [l2l . [l3| . In this setup, the gain medium 
can be thought of as a stream of suitably prepared atoms. 

From a conceptual point of view, a much simpler sys- 
tem relates to a single atom laser, where the gain medium 
is replaced by a single trapped atom. Such a laser has 
recently been experimentally demonstrated by Kimble's 
group [l6j , where a single atom interacts with a sin- 
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gle cavity mode. Thus the interesting question arises, 
whether a two-mode generalization of the single-atom 
laser also enables one to generate entanglement in macro- 
scopic light. 

Therefore, here we consider a single atom that interacts 
with two quantized modes of a doubly resonant cavity 
via two lasing transitions. In our model, the atomic level 
scheme is based on the single-atom laser experiment per- 
formed by Kimble's group [13], where dipole transitions 
between four hyperfine levels of atomic caesium were con- 
sidered. In contrast to their experiment, we do not work 
in the strong coupling regime since we are interested in 
the generation of large photon numbers. We show that, 
under certain realizable conditions, a two- mode single- 
atom laser can serve as a source of macroscopic entangled 
light. Macroscopic entanglement can be achieved over a 
wide range of control parameters and initial states of the 
cavity field. 

An important technical question in the generation 
of continuous variable entanglement in quantum opti- 
cal systems is the way such entanglement can be mea- 
sured experimentally. This is a hotly discussed subject 
in recent years. Several inequalities involving the cor- 
relation of the field operators have been derived that 
are based on the separability condition of the field 
modes 0, El, El S3 El S3, H, HI • A violation of these 
inequalities provides an evidence of entanglement. These 
inequalities can, in general, provide only a sufficient con- 
dition for entanglement and only, in some very specific 
instances, lead to sufficient and necessary conditions for 
entanglement. In this paper we use the inequality based 
on quadrature measurement of the field variables for the 
test of entanglement. 



II. MASTER EQUATION FOR THE DENSITY 
OPERATOR OF THE CAVITY MODES 



We consider a single four-level atom trapped in a dou- 
bly resonant cavity (see Fig.[T]). The atom interacts with 
two (nondegenerate) cavity modes and two classical laser 
fields. The intensities and frequencies of the two laser 
fields can be adjusted independently. The aim of this 
section is to derive an equation of motion for the reduced 
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FIG. 1: A single four-level atom is trapped in a doubly res- 
onant cavity and interacts with two cavity modes and two 
classical laser fields. The inset shows the atomic level scheme. 
The laser field with frequency U13 and Rabi frequency O3 
couples to the \a) <-> \d) transition, and the cavity mode 
with frequency v\ and coupling constant g\ interacts with 
the \a) <-> |c) transition. A a is the detuning of the fields ^3 
and gi with state \a). The laser field with frequency 014 and 
Rabi frequency Q4 drives the \b) <-> |c) transition, and the 
second cavity mode with frequency V2 and coupling constant 
#2 interacts with the \b) <-> \d) transition. A;, is the detuning 
of the fields SI4 and 52 with state \b). Spontaneous emission 
is denoted by dashed arrows, and the parameters Ti are the 
decay rates of the various transitions. 

density operator g-p of the two cavity modes. 

We begin with a detailed description of the system 
shown in Fig. [TJ The first cavity mode with frequency 
v\ couples to the atomic transition |a) «-> |c), and the 
second mode with frequency v 2 interacts with the atom 
on the 1 6) <-> \d) transition. In rotating- wave approxima- 
tion (RWA), the interaction of the atom with the cavity 
modes is described by the Hamiltonian 

H c = %iOi|o)(c| +hg 2 a,2\b)(d\ +H.c. (1) 

Here aj (aj) is the annihilation (creation) operator of 
the cavity mode with frequency Vj and coupling constant 
9j (j € {1,2}). The detuning of the first cavity mode 
with the \a) <-> |c) transition is denoted by A1, and A 2 
is the detuning of the second mode with the \b) «-> \d) 



transition, 

Ai = V\ - UJac , A 2 = V 2 - LObd ■ (2) 

The resonance frequencies on the \a) «-> \c) and \b) <-> \d) 
transitions have been labeled by oj ac and lom, respec- 
tively. In addition, the atom interacts with two classical 
laser fields. The first laser field with frequency CU3 and 
Rabi frequency O3 couples to the \a) <-* \d) transition, 
and the second field with frequency W4 and Rabi fre- 
quency $14 coherently drives the |6) «-> |c) transition. In 
rotating-wave approximation, the atom-laser interaction 
reads 

H L = -M7 3 |a)(d|e-^ 3t - hn 4 \b) (c\ e - luJit + H.c. . (3) 

Note that the Rabi frequencies ^3 = 1 exp(i</>3) and 
f^4 = ISI4I exp(i(/)4) are complex numbers, and <f>^ and 
(j>4 are determined by the phase of the laser fields. The 
detuning of the laser fields with the corresponding atomic 
transitions are 

A 3 = u 3 - u ad , A 4 = uj 4 - u bc , (4) 

where uj a d and io\,c are the resonance frequencies on the 
I a) «-» \d) and |6) <-> \c) transitions, respectively. 
The free time evolution of the cavity modes is governed 

by 

i/n = hvia\cii + hv2i\a,2 , (5) 

and Ha is the free Hamiltonian of the atomic degrees of 
freedom, 

H A = Hu a \a){a\ + tw b \b){b\ + Hlj c \c){c\ +fhu d \d){d\ . (6) 

With these definitions, we arrive at the master equation 
for the combined system of the atomic degrees of freedom 
and the two cavity modes, 

g=-^[H R + H A + H L + H c ,g}+C y g. (7) 

The last term in Eq. ([7|) accounts for spontaneous emis- 
sion and is given by 

1 4 

L -<e = ~ 2 E r * ( s t s ^e + esfsr ~ ^- g s+) , (8) 

where the atomic transition operators are defined as 

S+ = \a)(d\, S+ = \a){c\, 

S+ = \b)(c\, S+ = \b)(d\, Sr = (S+)i. (9) 

The parameters Ti are the decay rates of the various 
atomic transitions (see Fig. Q]). 

In a next step, we derive from Eq. |(7J) the master equa- 
tion for the density operator gp of the cavity modes, 

g F = Tr A g = g aa + Qbb + Qcc + Qdd , (10) 
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and q vv denotes (i/|p|z/). To this end, we apply a unitary 
transformation W — W R ® W\ to Eq. , where W R = 
exp[iHnt/h] acts only on the cavity modes, and 

W A =cxp[i{H A + fiA 3 |o)(o| + hA 4 \b)(b\)t/h] (11) 

acts only on the atomic degrees of freedom. As indicated 
in Fig. [TJ we assume that the condition of two-photon 
resonance is fulfilled, i.e. 

A a = Ai = A 3 , A b = A 2 = A 4 . (12) 

The density operator in the new frame is denoted by 
g = WgW^ and obeys the equation of motion 

g~ = -^[H + H c ,Q}+£ ig -, (13) 

where 

H Q = -hA a \a)(a\~hA b \b)(b\ 

- h(n 3 \a){d\ + 4 |&}(c| +H.c.) . (14) 

The two-photon condition Eq. (|12p ensures that the 
Hamiltonian Hq + He in Eq. (|13[) is time-independent. 
The master equation for the transformed density opera- 
tor gp of the cavity modes is obtained if we trace over 
the atomic degrees of freedom in Eq. (Ti"3"|) . 

qf = -ig\\a\, Qac] - ig2[a\, g~bd] + H.c. . (15) 

In order to eliminate the coherences g ac and Qbd from 
Eq. (fl5|) , we apply the standard methods of laser theory 
(see, e.g., Chapter 14 in 25]). We restrict the analysis to 



the linear theory and solve Eq. (fT3f to first order in the 
coupling constants g\ and g 2 . To this end, we expand the 
density operator gp in Eq. (fl"3|) as g = go + gc and retain 
only terms up to first order with respect to Hq. This 
procedure yields two uncoupled equations for go and gc, 

go =C go , (16) 
gc =C Q g c - z[ H c, Qo] , (17) 

and the superoperator £ is defined as 

Co{-) = ~[H ,-]+C 7 (-). (18) 

Here the centered dot denotes the position of the argu- 
ment of Co- The zeroth-order equation (fT6|) describes the 
interaction of the atom with the classical laser fields to 
all orders, and Eq. (fPfj) is the first-order equation. The 
steady state solution for g ac and g~bd can be obtained if 
the steady-state solution for go is plugged in Eq. (fl7|) . 
We find 

igiQac = auaiQF + ai 2 a\g F + /3ii^f«i + /?i20Fa 2 , 

ig2Qbd = a 22 a 2 g F + a 21 a\g F + (3 22 gFa 2 + foiQFat , 

(19) 

and the coefficients otij and are defined in Ap- 
pendix [£] Next we substitute Eq. (JT5J) in Eq. (JTSJ) to 
obtain the equation of motion for gp. Finally, we trans- 
form gp back with respect to Wr and obtain the equation 
of motion for the density operator gp of the cavity modes, 



OF = - ivi[a[ai, g F ] - iv 2 [a\a 2 , gp\ 



(20) 



a\\a\aigp + a^gpa^ai — (an + ajjjaigpaj — 

a 22 a 2 a 2 0F + a 22 £>F<4 a 2 ~ ( a 22 + Oi 22 )a 2 gpa 2 ~ 

(ai 2 + a 2 i)a\a\gp - (/3i 2 + P 2 i)gFa\a\ - (a 2 i 
(a* 12 + a 21 )g F axa 2 - (P$ 2 + Z^iWa^F - (a 21 
Ki (a\aigp + gpa\ai — 2aigpa\^j — k 2 {^a\a 2 gp ■ 



In the last line of Eq. (|2U)l , we included the damping of 
the cavity field. The damping constants of the cavity 
modes are denoted by K\ and n 2 , respectively. 

In the master equation (f2T)]) . the two classical laser 
fields are taken into account to all orders in the Rabi 



Plia\a\gp - (3ngpaia\ + (fin + ^n) a i^F a i 
I3 22 a 2 a\gp - f3 22 gFa 2 al + ((3 22 + f3 22 )alg F a 2 

- /3i 2 )a\gpal - (a 12 - I3 2 i)a\gpa\ exp[— i(y\ + v 2 )i\ 

- t3*i 2 )a 2 gpai - (a* 12 - l3 21 )aig F a 2 exp[£(i/i + v 2 )t] 

f gpa\a2 - 2a 2 gpa 2 
I 

frequencies ^3 and O4. On the contrary, the two quan- 
tum fields inside the cavity are only treated to second 
order in the coupling constants g\ and g 2 . This approxi- 
mation means that we ignore saturation effects and oper- 
ate in the regime of linear amplification. It is justified if 
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the Rabi frequencies associated with the quantum fields 
are small as compared to other system parameters which 
dominate the time evolution. 



III. ENTANGLEMENT OF THE CAVITY FIELD 



In this Section we show that the system depicted in 
Fig. [1] can serve as a source of macroscopic entangled 
light. We employ the sufficient inseparability criterion 
derived in [l8| to provide evidence for the entanglement 
of the two field modes. 

By definition, the quantum state qf of the cavity field 
is said to be entangled if and only if it is nonseparable, 
and g-p is separable if and only if it can be written as 



Qf 



E(i) 
Pie) 

3 



Q 



(2) 



(21) 



Here gp and g^ are normalized states of the modes 1 
and 2, respectively, and the parameters jp, > comply 
with Pi = 1- The criterion derived in [181 ] states that 
the system is in an entangled quantum state if the to- 
tal variance of two Einstein-Podolsky- Rosen (EPR) type 
operators u and v of the two modes satisfy the inequality 



((Au) 2 + (Av) 2 ) <2. 



where 



X\ + x 2 



v = pi - p 2 . 



(22) 



(23) 



Here ik and pk are local operators which correspond to 
mode k with frequency . They must obey the commu- 
tation relation 



[x k ,pi] = iS k i , 



(24) 



but are otherwise arbitrary. For the physical system con- 
sidered here, it turns out that the following quadrature 
operators 

x k = (bk + b\)/V2 and Pk = {b k - b\)/(V2i) (25) 
are the best choice, where 
bk(t) = etfc exp[ii/fci] and b\,(t) = a\. exp[— ivkt] ■ (26) 



With the help of Eqs. (|23)) and (|25|) . we express the total 
variance of the operators it and v in terms of the operators 
bk and fot, 



( (Aw) 2 + (At)) 2 > = 2 [l + <6+6i> + (b\b 2 ) + (hh) 

+{b\b\) - (h){b\) (b 2 )(bl) (&x)(6 2 ) - (b\)(bl) . 

(27) 

In Appendix [B] we outline the calculation of the mean 
values that enter Eq. ([2T]) . 
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FIG. 2: (a) Time evolution of {(Aft) 2 + (AC) 2 ). The mean 
value of the total number of photons {N) is shown in (b) on a 
logarithmic scale. At t = 0, the cavity field is assumed to be 
in the vacuum state. The dashed curves were obtained with 
the density operator of the parametric oscillator in Eq. (|30[) . 
and the solid curves correspond to the full density operator 
in Eq. I] 20 [I. The parameters are gi — g2 — g, IO3 1 = 25g, 

|fi 4 | = 2 g , ih = r 2 = r 3 = r 4 = 5<?, A a = 0, A b = 40<?, 



K2 = 10 g and 4>i + 4>i = tt/2. 



Next we classify several parameter regimes for which 
the inequality ([22]) is fulfilled. In a first step, we consider 
the case where the Rabi frequency and the detuning 
At are much larger than the parameters |A a |, | 1 , Fj 
(ie {1,2,3,4}), i.e. 



|fi 3 |, |A 6 | » IAJ, |fi 4 |, Ti. 



(28) 



If these conditions are fulfilled, the parameters and 
in Eqs. (|AT] ) -(|A8] ) of Appendix O reduce to 

an PS 0, a 22 ~ 0, n ~ 0, f3 2 2 ~ , 

a 2 i ps 0, fa PS 0, ai2 ps (i 2 i ps ~iaexp[z(0 3 + cj>i)t] , 

in 4 | 



a = 5i32 ■ 



fi a |A h ' 



(29) 
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In these equations, <j>% and (f>4 are the phases of the classi- 
cal laser fields with Rabi frequencies S! 3 = |fi 3 | exp(itp 3 ) 
and | 1 exp(i04), respectively (see Sec. Ql]). If the ap- 
proximate parameters in Eq. (|29p are plugged in Eq. (j2"0")l . 
we obtain the equation of motion for the density operator 
qf of the cavity modes in the limit (|28|) , 



T. 



z ^i[aia 1; p F j - i v 2 W 2 a2, Qf] + i[Hp 
k (a\aiQF + QFa\ai — 2aiQFa\ 

gFd\a2 — 2q2Qfo-2 



a\a,2QF 



(30) 



where 
Hp =aa\a) 2 eyi]>[i 



3 + 4 )f] exp[—i(v-L + V2)t] 

aaiCL2 exp[— i(4>3 + 4>i)t] exp[i(i^i + V2)t\ 



(31) 



Here we assumed for the sake of simplicity that the de- 
cay rates of the cavity modes are equal, k\ = k-2 = n. 
We identify Eq. ([30]) as the master equation for a nonde- 
generate parametric oscillator in the parametric approx- 
imation [251 ] . Note that this parametric limit was also 
obtained in the case of a two-mode correlated sponta- 
neous emission laser discussed in [l2j|. Next we evaluate 
the total variance of the operators u and v in Eq. |(77|) 
and the mean number of photons (N) = (a[ax + a\a 2 ) = 
{b\b\ + 62^2) with the approximate density operator qf in 
Eq. p0|) . If the sum of the laser phases obeys <p3 + <j>4 — 
7r/2, we obtain [l2j 



((AuY + (AvY)(t)= ((AuY + (Av) 2 )(0) 



2 k 



2 k 



e -2( a+K )t + (32) 

a + K 



(N)(t) 



(N)(0) 



cosh(2ai)e 



+ (a 1 a 2 + a\al)(0) 



-iKt 



sinh(2o:i)e 



-Int 



or 



(33) 



It follows from Eq. (j3"2"|) that the entanglement criterion 
in Eq. ([22]) is satisfied for any initial state of the cavity 
field if (a + n)t 3> 1 and a > [l2j . The time evolution 
of the total variance of the operators u and v is shown in 
Fig. Ufa). The dashed curve shows ((Aw) 2 + (Av) 2 ) ac- 
cording to Eq. ([32]) , and the solid line corresponds to the 
general case where the mean values in Eq. (|27j) are eval- 
uated with the full density operator qf in Eq. (|2"0"]) . The 
cavity modes are assumed to be in the vacuum state ini- 
tially, and the parameters comply with condition (|28p . 
It follows from Fig. [2] that the approximate result in 
Eq. (|3"2"|) is only in good agreement with the exact so- 
lution if gt < 300. While the light field remains in an 
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FIG. 3: (Color online) (a) Time evolution of ((Ai) 2 + (Av) 2 ). 
The mean value of the total number of photons (N) is shown 
in (b) on a logarithmic scale. At t = 0, the cavity field is 
assumed to be in the vacuum state, and we set Fi = T2 = 
r 3 = T 4 = 5g, gi = gi = g, Ki = k 2 = W 3 g and <j) 3 + 
4>4 — 7r/2. The parameters for the curves labeled with I are 



l«3| : 
have 



= 25 ff , |n 4 | 
n 3 \ = 15g, 



= 9.8<?, A a = 
fi 4 = Gg, A a 



0, A;, = 43g, and for II we 
= 0, A b = 32.5g. 



entangled state in the parametric case, the exact solu- 
tion demonstrates that the entanglement of the cavity 
field exists only for a finite period of time. 

Next we discuss the time evolution of the mean num- 
ber of photons (N). According to Eq. (fBU)) . (N) grows 
exponentially with time for any initial state of the cavity 
field, provided that (a — n)t 1 and a > n [12]. The 
time evolution of (N) is shown in Fig.[2]Jb) on a logarith- 
mic scale. In contrast to ((An) + (Av) ), the result for 
(N) in the parametric approximation (dashed line) is in 
good agreement with the exact solution (solid line) even 
for gt ^> 300. Moreover, Fig. [DJb) shows that the mean 
number of photons grows exponentially if the scaled time 
gt is sufficiently large. 

According to Fig. ® the entangled state of the cavity 
field contains up to (N) « 110 photons on average. It 
follows that the single-atom laser depicted in Fig. [1] can 



6 



give rise to an entangled quantum state of the two cavity 
modes if the parameters are in agreement with condi- 
tion ([28]) . If this condition holds, level \b) is almost not 
excited due to the large detuning At,, and states |c) and 
\d) are coupled via a two-photon process. In contrast, the 
transitions \d) <-> \a) and |c) <-> \a) are driven resonantly. 
In this situation, the structure of the Hamiltonian Hp in 
Eq. (|3ip implies that the system can only emit photons 
into the cavity fields in pairs, where one photon is emit- 
ted in mode 1 and the other photon in mode 2. If the 
cavity field is initially in the vacuum state 1 0, 0) , it will 
evolve under the influence of Hj> into the entangled state 

a|0,0) + 6|l,l) + c|2,2) + ... , (34) 

where a, b and c are complex coefficients. If the compli- 
cated master equation ([20)l can be reduced under certain 
conditions to the parametric equation (I30j) , it is thus clear 
that a macroscopic entangled state is generated. 

Due to the symmetry in the atomic level scheme, it is 
possible to reverse the role of the transitions \d) *-* \a) «-> 
|c) and |c) <-> \b) <-> \d). In this case, the detuning A a 
is large and the transitions \d) <-> \b) and |c) <-> \b) are 
driven resonantly. Condition (|28|) then has to be replaced 
by 

|fi 4 |, |A„| » |A b |, |fi 3 |, Ti, (35) 

and the only nonvanishing coefficients in Eq. (|20j) are now 
determined by an ~ (3\2 ~ — ia' exp[i(</>3 + <p4)t], where 
a' = gig2|^3|/(|^4|A a ). It follows that the results in 
Eqs. ([50]). p2|) and (J33J) are also valid if condition (f35|) 
holds, provided that a is replaced by a' . 

We now demonstrate that it can be advantageous to 
consider parameters which do not comply with condi- 
tions ([28]) or (f35|) . Since the approximate results in 
Eqs. (|32p and (|3"3"|) do not apply in this case, we eval- 
uate the mean values ((Am) 2 + (Aw) 2 ) and (N) only with 
the exact density operator qf in Eq. (|20|) . The time evo- 
lution of ((Aii) + (Aw) ) is shown in Fig. 02a) for two 
sets of parameters. As compared to the parameters cho- 
sen for Fig. ([2]), the magnitude of the Rabi frequency f^4 
has been increased such that |Oa| is still larger, but not 
much larger than \Q.4\. It follows from Fig. [3]^ a) that the 
entanglement criterion in Eq. ([2"2"[) is fulfilled for shorter 
times as compared to the solid line in Fig. [Ha). On the 
other hand, Fig. El^b) shows that the mean number of 
photons can be much larger as compared to Fig. [U^b). 
For curve I of Fig. [31[a), the maximum mean number of 
photons for which the entanglement criterion ([22]) is still 
fulfilled is (N) ps 10.2 x 10 4 . The same number for the 
parameters of curve II reads (N) w 6100. As compared 
to Fig. [21 the maximum mean number of photons can be 
enhanced by several orders of magnitude. 

Finally, we consider the case where the quantum 
state of the cavity field is initially the coherent state 
[100, -100). The time evolution of ((Aii) 2 + (Aw) 2 ) and 
(N) is shown in Fig. Q] for two sets of parameters. All 
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FIG. 4: (Color online) (a) Time evolution of ({Aii) 2 + (Avf). 
The mean value of the total number of photons (N) is shown 
in (b) on a logarithmic scale. At t = 0, the cavity field is 
assumed to be in the coherent state 1 100, —100), and we set 
Ti = T 2 = T 3 = T 4 = 2g, 51 = g 2 = g, «i = n 2 = 10~ 2 g and 
4>s + "^4 = tt/2. The parameters for the curves labeled with I 
are | $^3 1 = 10<?, |fi4| = 5g, A a ~ 0, Af, = 15g, and for II we 
have |n 3 | = Wg, |^4| = 2g, A a = 0, A b = I5g. 



mean values were evaluated with the exact density oper- 
ator in Eq. ([2"0"[) . For curve I, the magnitude of the Rabi 
frequency ^4 is larger as compared to curve II. All other 
parameters are the same for curve I and II. It can be seen 
from Fig. 0Ja) that the entanglement criterion is fulfilled 
for shorter times if |fi4| is increased. In contrast, the 
mean number of photons can be greatly enhanced if the 
value of 1^4 1 is increased, as can be seen from Fig. Sfb) . 
Similar conclusions can be drawn from Fig. (3[ where the 
initial state of the cavity field is the vacuum. The com- 
parison of Figs. [3] and [4] shows that the mean number 
of photons can be much larger than in Fig. [3] if the cav- 
ity field is initially prepared in a coherent state. Due to 
the large mean number of photons in the cavity modes, 
the system may leave the regime of linear amplification 
such that saturation effects modify the curve progression 
in Fig. Q] These effects are described by terms that go 
beyond the second-order expansion of the atom-cavity 
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coupling and are neglected here. According to the lin- 
ear theory, the maximum mean number of photons for 
which the entanglement criterion (|22|) is still fulfilled is 
(N) « 6.5 x 10 5 in the case of curve II of Fig. 2J and in 
the case of curve I the entangled cavity field contains up 
to (N) m 5.4 x 10 7 photons. 



IV. CONCLUSION 

We have shown that a two-mode single-atom laser can 
serve as a source of macroscopic entangled light. We iden- 
tified two parameter regimes for which the quantum state 
of the cavity field is entangled for a long period of time. 



For these parameters, the master equation for the density 
operator of the two cavity modes can be approximately 
reduced to the master equation for a nondegenerate para- 
metric oscillator in the parametric approximation. 

The mean number of photons in the cavity field can be 
strongly increased if parameters beyond the parametric 
limit are chosen. This enhancement of the mean pho- 
ton numbers is accompanied by a shortening of the time 
slice for which the entanglement criterion is fulfilled. As 
the initial state of the cavity field, we chose either the 
vacuum or a coherent state. We demonstrated that the 
mean number of photons of the entangled cavity field can 
increase by several orders of magnitude if a coherent state 
instead of the vacuum is chosen as an initial state. 



APPENDIX A: COEFFICIENTS 

Here we give the explicit definitions of the coefficients otij and /Jy which enter the master equation (|20p for the 
density operator g-p of the two cavity modes 



a n =2g 2 T 2 |ft 3 | 2 |ft 4 | 2 [4 (P 2 * + AiA b ) |ft 4 | 2 + P 4 * (a \Sl 3 \ 2 + P 1 (P 1 + P 2 *))] /(P 3 P 4 ) , 

fa = - 2g 2 1 T i |f» 3 | 2 |ft 4 | 2 [4Pi |0 4 | 2 + P 2 * (4 |f! 3 | 2 + P 1 (P 1 + P 2 *))] /(P 3 P 4 ) , 

au = - 2g l92 T 2 n 3 n 4 \n 3 \ 2 [4^ \n 3 \ 2 + P x 2 (P + P 2 *) - 4 |ft 4 | 2 (2P a + P*)] /(P 3 P 4 ) , 

p 12 = - 2g l92 T i Q a n i |ft 4 | 2 [(P + P 2 *) |P 2 | 2 + 4 |f! 4 | 2 P 2 + 4 \Q 3 \ 2 (P - 4zA a )] /(P 3 P 4 ) , 

a 22 =2. 9 2 r 4 |r! 3 | 2 |r! 4 | 2 [4(P*+4iA Q )|f! 3 | 2 + P 2 * (4 |fi 4 | 2 + P 2 (P 2 + P*))] /(P3P5) , 

(3 22 = - 2g\Y 2 |ft 3 | 2 |r» 4 | 2 [4P 2 |r» 3 | 2 + P* (4 |r» 4 | 2 + P 2 (P 2 + P*))] /(P3P5) , 

a 2 i = - 2g l92 T i n 3 Q i |fi 4 | 2 [4P 2 |fi 4 | 2 + P 2 (P 2 + P*) - 4 \<TL 3 \ 2 (2P 2 + P*)] /(P3PO , 

p 21 = - 2g 1 g 2 r 2 n 3 n i \n 3 \ 2 Up 2 + p*) |p| 2 + 4 |f> 3 | 2 p + 4 |n 4 | 2 (p 2 - 4*a 6 )1 /(p 3 p 5 ) . 



The parameters Pi, P 2 , P 3 , P 4 and P5 in Eqs. (|Al|) - (|A8|) are defined as 



Pi 


=r 3 


+ r 4 


+ 2iA b , 










Pi 


=r 4 


+ r 2 


+ 2iA a , 










P3 


=r 2 


P1I 2 


n 3 \ 2 + r 4 




+ 8 


^ 3 | 2 l«4| 2 (r 2 4 


r 4 ) 


Pi 


-< 




H^l 2 ) 2 


+ Pi (P + 




\n 3 1 2 + p 2 *(Pi 


+ p; 


P 5 


-< 


l^sl 2 


-|^ 4 | 2 ) 2 


+ Pi (P 2 ■+ 


P* 


)|n 3 | 2 + p 2 (p 2 


+ Pi 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
(A8) 

(A9) 
(A10) 
(All) 

(A12) 
(A13) 



APPENDIX B: CALCULATION OF THE MEAN VALUES 



In the following, we outline the calculation of the mean values that enter the total variance of the operators u and 
v in Eq. (|2~?| . We begin with the mean values of the quadrature operators defined in Eq. (f2"S"]) with respect to the 
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density operator g-p of the two cavity modes. With the help of Eq. (|20|) , we derive the following system of differential 
equations for the mean values (bi) and (b 2 ), 



(bi) 
(bl) 



Cli + K\ C\ 2 



(bl 



21 



ch + «2 / V (4) 



and Cij = on 



Pij . The solution to this set of coupled equations is 



(bi) =6 



(6 



t\ —„w 2 t 



cosh(iui£)(&i) + ((&i)o (C 22 - Cn - «i + k 2 ) - 2(4)oCi 2 ) sinh(wit) 
cosh(u; 1 t)(6 2 )o + — ((fe 2 )o (C u - C* 2 + Ki - k 2 ) - 2(6 1 ) C* 1 ) sinh(wit) 



(Bl) 

(B2) 
(B3) 

., ^C X2 C* 2X + (C n - C* 22 + ki - k 2 ) 2 , w 2 = -|(C 11 + C* 2 + Kl + K2 ), (B4) 

and ( • )o = (-)(i = 0) denotes the initial mean value at t = 0. Note that the mean values (&{) and (62) can be 
obtained from (£>i) and (bl) by complex conjugation, i.e. (b\) — (b\)* and (b 2 ) — (b 2 )* . 

The remaining mean values in Eq. ([27]) involve products of the operators bi and bj . With the aid of Eq. (f2"0"l) , we 
obtain the following set of differential equations, 

d t R = MR + I , 



where 



Wl 



(B5) 



and 



where R = (^(b[bi), (b\b 2 ), (hh), (b\b\ 

( Dn 
D 22 



M = - 



Cf 



C12 \ 



C21 C\ 



C 2 i C21 
D 12 

D* 12 J 



( fax + 0ti \ 
P22 + Ph 
ol\2 + a 2 i 

\ a*2 + "21 / 



(B6) 



V C* 21 C* 12 

The elements of the matrix M are defined as 

Cij = aij + fa , D u = a u + a*u + Pu + /?* + 2 Ki , and D 12 = Cn + C 22 + Ki + k 2 ■ (B7) 

The differential equation Eq. (|B5|) can be solved numerically without difficulties. An analytical solution can be 
obtained, for example, by means of the Laplace transform method which yields the following results for the components 
Ri of the vector R, 

4 

Ri = J2 [Res(/ t , A fc ) + Res( ffl , A*)] e Xkt + Res( ffl , 0) . (B8) 

k=l 

In this equation, expressions of the type Res(/i, z) denote the residue of the function h evaluated at z, and the functions 
/ = (/1, / 2 , /3, fi) and g = (g x , g 2 , g 3 , g 4 ) are determined by 

/(*) = [s U ~ M] ~* R a and g(s) = [s 1 4 - M] _1 (I/s) , (B9) 

respectively. Here I4 denotes the 4x4 identity matrix, and the vector Rq is the initial value of R at t = 0, 
Rq = \{b\bi)a, (b 2 b 2 )o, (b±b 2 )o, (b\b\)o^. Finally, the parameters Afc are the four (complex) eigenvalues of the matrix 
M which is defined in Eq. (|B6[) . The eigenvalues Afc can be obtained as the roots of the following equation, 
s 4 + [Dn + D 22 + D 12 + D$ 2 ] s 3 + [\D 12 \ 2 - 2C 21 C* 12 ~ 2C 12 C* 21 + D 1X D 22 + (Dn + D 22 )(D 12 + D* 12 )] s 2 

+ [(Dn + D 22 )\D 12 \ 2 ~ (C 2X Cl 2 + C 12 C* 21 )(Dn + D 22 + D 12 + D* 12 ) + D n D 22 (D 12 + D* 12 )] s 

+ C 2 2l (C* l2 ) 2 - (2C 12 C* 2l + DnD 12 + D 22 D* 12 )C 21 C* 12 + (C 12 C* 21 - D 22 D 12 )(C 12 C* 2l - DnD* 12 ) = . (BIO) 
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